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Abstract: Recently, the concept of tail dependence has been discussed in financial applications 
related to market or credit risk. The multivariate extreme value theory is a proper tool to measure 
and model dependence, for example, of large loss events. A common measure of tail dependence 
is given by the so-called tail-dependence coefficient. We present a simple estimator of this latter 
that avoids the drawbacks of the estimation procedure that has been used so far. We prove strong 
consistency and asymptotic normality and analyze the finite sample behavior through simulation. 
We illustrate with an application to financial data. 
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1 Introduction 

■ Modern risk management is highly interested in assessing the amount of extremal depen- 

dence, a growing phenomena in recent time periods of volatile and bear markets. Cor- 
' relation itself is not enough to describe a tail dependence structure and often results in 

misleading interpretations (see e.g. Embrechts et al. [3], 2002 for examples). 



Multivariate Extreme value theory (EVT) is the natural tool to deal with tail depen- 
dence. The so-called tail- dependence coefficient (TDC) has become a popular measure in 
OO . risk management. It is usually denoted A and measures the probability of occurring ex- 

' treme values for one random variable (r.v.) given that another assumes an extreme value 

• • too. More precisely, 

A = limP(Fx(X) > l-t|Fy(y) > l-t), (1) 

where Fx and Fy are the distribution functions (d.f.'s) of r.v.'s X and Y, respectively. 
^ . The TDC characterizes the dependence in the tail of a random pair {X,Y), in the sense 

^ I that, A > corresponds to tail dependence whose degree is measured by the value of A, 

and A = means tail independence. Due to the emergent importance of this issue, it is 
not surprising that the implementation of tail dependence measures and respective estima- 
tion have attracted the attention of investigators. Sibuya (|17|. 1960), Tiago de Oliveira 
([IB], 1962-63), Ledford and Tawn ([ISKH, 1996, 1997), Joe ([I^, 1997), Coles et al. ([2], 
1999), Embrechts et al. (|3|, 2003), Frahm et al. ([8], 2005), Schmidt and Stadtmiiller 
(|16|. 2006), Ferreira and Ferreira ([6], 2011), are some references on this topic. Curiously, 
the first tail dependence concept appearing in literature concerns the TDC, as far back in 
the sixties with Sibuya (|17|. 1960). More precisely, for Normal distributed random pairs, 
Sibuya (|17|. 1960) shows that no matter how high we choose the correlation, if we go far 
enough into the tail, extreme events appear to occur independently in each margin. 
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Here we shall present an estimator for the TDC based on a new procedure that avoids 
the main drawback of existing ones. Strong consistency and asymptotic normality are 
stated and the finite sample behavior is illustrated through simulation. An application to 
financial data is presented at the end. 

2 EVT and tail dependence 

The main objective of an extreme value analysis is to estimate the probability of events 
that are more extreme than any that have already been observed. The tools within EVT, 
in particular the use of extremal models, enables extrapolations of this type. The cen- 
tral result in univariate Extreme Value Theory (EVT) states that, for an i.i.d. sequence, 
{Xn}„>i, having common distribution function (d.f.) F, if there are real constants a„ > 
and bn such that, 

P(max(Xi, ...,Xn) < ttnX + bn) = F'^ianX + bn) >n^co G{x) , (2) 

for some non degenerate function G^, then it must be a Generalized Extreme Value function 
(GEV), 

G{x) = exp(-(l + -fx)-^/^), 1 + 72; > 0, 7 G M, 

(for 7 = 0, G{x) = exp(— e~^')) and we say that F belongs to the max-domain of attraction 
of G, in short, F £ T){G). The parameter 7, known as the tail index, is a shape parameter 
as it determines the tail behavior of being so a crucial issue in EVT. More precisely, 
if 7 > we are in the domain of attraction Frechet corresponding to a heavy tail, 7 < 
indicates the WeibuU domain of attraction of light tails and 7 = means a Gumbel domain 
of attraction and an exponential tail. 

Models for dependent or non-identical distributed r.v.'s have also been developed (see, 
for instance, Leadbetter |12| 1983 and Mejzler |15| 1956, respectively). 

In a multivariate framework, an extension of the univariate limiting result in ([2]) is 
considered and the multivariate maximum corresponds to the vector of component-wise 
maxima. Let {Xn = (X„^i, X„^rf)}„>i be an i.i.d. sequence of d-dimensional random 
vectors with common d.f. F. If there are real vectors of constants bn and positive an such 
that, 

P(max(Xi,...,Xn) < anX + bn) = F"(anX + bn) ^„^oo G'(x), (3) 

for some non degenerate function G, then it must be a multivariate extreme value distri- 
bution (MEV), given by 

G(x) = exp(-/(- log Gi(xi), - log Gd{xd))), 

for some d-variate function /, where Gj, j = 1, is the marginal d.f. of G. We also say 
that F belongs to the max-domain of attraction of G, in short, F G T){G). The function I 
in is called stable tail dependence function. 

Let Fj be the marginal d.f. of F. Since a sequence of random vectors can only converge 
in distribution if the corresponding marginal sequences do, we have, for j = 1, d, 

Pj {^n,jXj -\- bnj) ^n— >oo Gj{Xj) 
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Hence Gj is a GEV and Fj is in its domain of attraction. 

In order to study the dependence structure of a MEV, it is convenient to standard- 
ize the margins so that they are ah the same. A particular useful choice is the unit 
Frechet, exp(— (observe that unit Frechet marginals can be obtained by transforma- 
tion — 1/log Fj{Xj) for j S / C {1, d}), and the stable tail dependence function / in (j4]) 
becomes 

/(v) = -logG(^;r\...,t;^i), vG [0,oo). (4) 

Thus I satisfies an important homogeneity property (of order 1) and hence, it is easy to 
establish that, except for the special case of independence, all bivariate extreme value 
distributions (BE V) are tail dependent (A > 0). Furthermore, we have 

A = 2-/(1,1). (5) 

Examples of parametric BEV models 

• Logistic: l{vi,V2) = (f|/r + V2/ry , with Vj > and parameter < r < 1; complete 
dependence is obtained in the limit as r — )■ and independence when r = 1. 

• Asymmetric Logistic: l{vi,V2) = (1 — ti)vi + (1 — t2)v2 + {{tiviy/r + {t2V2)^/rY, 
with Vj > and parameters < r < 1 and < tj < 1, j = l,2; when ti = t2 = 1 
the asymmetric logistic model is equivalent to the logistic model; independence is 
obtained when either r = 1, ti = or t2 = 0. Complete dependence is obtained in 
the limit when ti = t2 = 1 and r approaches zero. 

• Hiisler-Reiss: l{vi,V2) = vi^{r~^ + log(t;i/t;2)) + V2^{r~^ + log(t;2/f i)), with 
parameter r > and where $ is the standard normal d.f.; complete dependence is 
obtained as r — )• oo and independence as r — t- 0. 

3 Estimation 

The (i-variate stable tail dependence function in (jH) can also be formulated as 

hm tp(F,{X,) > 1 - a V ... V Fa{X,) > 1 - ^) (6) 

since, by applying the unit Frechet marginals, we have successively, 

limt^^tp(Fi{Xi) > 1 - ^ V ... V FdiXd) >l-f) 
= limt^oo tptxi > J- V ... V Xrf > ^ 
= limt^oo tPi^Xi > J- V ... V Xrf > ^ 

= lim*^oo-logF*(;^,...,|^) 

Therefore, based on ([5]) and ([7]), the TDC in ([1]) can be stated like follows: 

A = 2 - lim tp(Fi{Xi) > 1 - - V ^2(^2) > 1 - -V (8) 

t^oo \ t tJ 
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Huang (1992 [9]), considered the estimator based on ([8]) by plugging-in the respective 
empirical counterparts, 



i=l 



where Fj is the empirical d.f. of Fj, j = 1,2. Concerning estimation accuracy, some 
modifications of this latter may be used, like replacing the denominator n by n + 1, i.e., 
considering 



^ 1 



n + _ 

fe=i 

(for a discussion on this topic see, for instance, Beirlant et al. [1] 2004). The consistency 
and asymptotic normality of estimator A^''^^ is derived under the condition that is an 
intermediate sequence, i.e., fe^ — >• oo and kn/n — t- 0, as n — >• oo. The choose of the value k 
in the sequence {kn} that allows the better trade-off between bias and variance is of major 
difficulty, since small values of k come along with a large variance whenever an increasing 
k results in a strong bias. Therefore, simulation studies have been carried out in order to 
find the best value of k that allows this compromise. The other estimators, either for the 
stable tail dependence function I or for the TDC, that have been considered in literature 
(for a survey, see for instance, respectively, Krajina |11| 2010 and Frahm et al. |8] 2005) 
are also based on asymptotic results with the same drawback of including an intermediate 
sequence, already referred above. 

The approach that is presented here avoid this problem. It is based on an estimation 
procedure for the stable tail dependence function only involving a sample mean. More 
precisely, by Proposition 3.1 in Ferreira and Ferreira (0, 2011), we have that 

, _ j^(Fi(Xi)V"WF2(X2))V"^ 

i[XuX2)- vF2(X2))V-2- ^'^^ 

Therefore, based on ([5]) and (|10p we propose estimator 



A = 3-(l-Fi(Xi) VF2(X2))-^ (11) 



where Fi{Xi) V ^2(^2) is the sample mean of Fi{Xi) V ^2(^2), i.e., 
Fi(Xi) V F2{X2) = -Y^ [Fi{xl''>) V F2(XW)] . 



n 

i=l 



Proposition 3.1. Estimator \ in is asymptotically normal whenever F has continu- 
ous marginals and continuous partial derivatives. Moreover it is strong consistent. 

Dem. By Fermanian et al. ([S], 2002, Theorem 6), we have that 
^^{J(Fi(xP),F2(X«))-i?(J(Fi(^P),F2(^?^)))}^ / Giu,,U2)dJiuuU2) 



in distribution in £°°{[0, 1]^), where the limiting process and G are centered Gaussian, and 
J{ui,U2) = max(t/i, ti2), (ui,it2) G [0,1]^. The asymptotic normality is now derived from 
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a general version of the Delta Method as considered in Schmidt and Stadtmiiller |16| (2006; 
Theorem 13). 

The strong consistency is straightforward from Proposition 3.7 in Ferreira and Ferreira 
(0,2011). □ 

Remark 3.1. Observe that, if the marginals Fj, j = 1, 2, in are known, the asymptotic 
normality of X is straightforward by the Central Limit Theorem and the usual Delta Method. 
More precisely, 

V^(A- A) ^ iV(0,(j2), (12) 

where 

2 ^ ;(i,i)(i+;(i.i))' 
^ (2+Ki,i)) 

Details of the calculations can be seen in Proposition 3.3 of Ferreira and Ferreira (JTj, 
2011). The strong consistency is also immediately derived from the sample mean. 

3.1 Simulations 

We consider 1000 independent copies of n = 50, 100, 500, 1000 i.i.d. pseudo-random vectors 
generated from three different BEV models considered in Section O logistic, asymmetric lo- 
gistic and Hiisler-Reiss. We estimate the TDC through our estimator (A). For comparison, 
we also compute estimator A^^^^ and, the required choice of k to balance the variance-bias 
problem is based on the procedure in Schmidt and Stadtmiiller (|16). 2006). The empirical 
bias and the root mean-squared error (rmse) for all implemented TDC estimations are 
derived and presented in Table [TJ Our estimator A clearly outperforms estimator A*^^-*. 

3.2 Application to financial data 

We shall see evidence of tail dependence in financial data. We consider the monthly 
maximum of the negative log-returns of Dow Jones and NASDAQ index for the time period 
1994-2004. The corresponding scatter plot and TDC estimate plot of A^^^ for various k 
(Figure [1]) show the presence of tail dependence and the order of magnitude of the tail- 
dependence coefficient. Moreover, the typical variance-bias problem for various threshold 
values k can be observed, too. In particular, a small k induces a large variance, whereas 
an increasing k generates a strong bias of the TDC estimate. The threshold choosing 
procedure of k used in Section [3] leads to a TDC estimate of A^^^ = 0.5556 and from our 
estimator we derive A = 0.5268. 
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Table 1: Sample absolute bias and RMSE of the non-parametric TDC estimators A and 
for BEV Logistic, Asymmetric Logistic and Hiisler-Reiss models. 
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A 
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Figure 1: Scatter plot of Dow Jones versus NASDAQ monthly maximum negative log- 
returns (n = 84 data points) and the corresponding TDC estimates A^^-* for various k/n. 
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